Dynamical Locking of the Chiral and the Deconfinement Phase Transition in QCD 
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We study the fixed-point structure of four-fermion interactions in two-flavor QCD with 7V C colors 
close to the finite-temperature phase boundary. In particular, we analyze how the fixed-point 
structure of four-fermion interactions is related to the confining dynamics in the gauge sector. We 
show that there exists indeed a mechanism which dynamically locks the chiral phase transition to 
the deconfinement phase transition. This mechanism allows us to determine a window for the values 
of physical observables in which the two phase transitions lie close to each other. 
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I. INTRODUCTION 

The relation of quark confinement and chiral symmetry 
breaking in quantum chromodynamics (QCD) is not yet 
fully understood. In particular at finite temperature and 
quark chemical potential the investigation of the QCD 
phase boundary represents one of the major research top- 
ics in theoretical physics and is also of great importance 
for a better understanding of heavy-ion collision experi- 
ments [l| . Theoretical studies in the limit of many colors 
in fact suggest that an understanding of the interrelation 
of the chiral and the deconfinement phase transition is 
required to comprehend the QCD phase structure close 
to a possible critical endpoint @ ■ 

While the confinement transition is driven by the gauge 
degrees of freedom, the chiral phase transition is trig- 
gered by strong quark self-interactions. In Nambu-Jona- 
Lasinio (NJL) models (and in the Polyakov-loop ex- 
tended version thereof) these quark interactions are con- 
sidered as parameters 0-E3] and used to fit a given set 
of low-energy observables. In full QCD, however, we ex- 
pect that the quark self-interactions are dynamically gen- 
erated and driven to criticality by the gauge degrees of 
freedom. This has indeed been confirmed by means of a 
renormalization group (RG) analysis of the influence of 
gluodynamics on the fixed-point structure of four-fermion 
interactions, see e. g. Refs. [l5l - [2ll ]: once the gauge cou- 
pling exceeds a critical value, the quark sector is driven 
to criticality without requiring any fine-tuning. This ob- 
servation already suggests that there might be a deeper 
relation between the chiral dynamics in the matter sec- 
tor and the confining dynamics in the gauge sector and 
serves as a motivation for the present study. 

The deconfinement phase transition has been studied 
in pure SU(N C ) gauge theories and in QCD with lat- 
tice simulations, see e. g. Refs. I22h33 |. as well as func- 
tional continuum methods [33M40j | . In so-called PNJL 
and Polyakov-loop extended quark-meson (QM) models 
a background field (Aq) is introduced to study some 
aspects of quark confinement and the associated phase 
transition. This background field can be related to the 
so-called Polyakov variable Lp4o], 



where j3 is the inverse temperature, N c is the num- 
ber of colors, g denotes the bare gauge coupling and 
V stands for path ordering. In fact, it has been shown 
that trp L[(Ao)] serves as an order parameter for quark 
confinement in Polyakov-Landau-DeWitt gauge [33, H3| 
where (Aq) is an element of the Cartan subalgebra 
and denotes the ground state of the associated order- 
parameter potential in the adjoint gauge algebra. 1 This 
potential can be computed, e. g., from the knowledge of 
gauge correlation functions. In Fig. [1] we show the re- 
sults for the order-parameter potential as obtained from 
a first-principles RG study [H, HI] . The order parame- 
ter tvp L[(A )} is related to the standard Polyakov loop 
(trp L[j4q]) y i a the Jensen inequality, 



tT F L[(A )] > (tr F L[A)]) 



(2) 



In PNJL/PQM model studies one of the underlying 
approximations is to set trpL[(^o)] = {^fL[Ao}). This 
opens up the possibility to incorporate results for the 
Polyakov loop (trpLpAo]) as obtained from lattice sim- 
ulations in these studies It is then found that 
the chiral and the deconfinement phase transition lie in- 
deed close to each other at small values of the quark 
chemical potential, as one would naively expect it to 
be the case in full QCD [HM M, M^M- However, 
the inclusion of the Polyakov loop in PNJL/PQM model 
studies by means of a Polyakov-loop potential is not 
unique and the QCD phase boundary at finite chemi- 
cal potential has indeed been found to be very sensitive 
to different parameterizations of this potential, see e. g. 
Ref. |4l|. On the other hand the data for the order pa- 
rameter L[(Ao)] is available. Therefore it seems natural 
to study at least some of the consequences of the approxi- 
mation tip L[(Aq)] = (trp L[Aq\) underlying these model 
studies. 

Recently, so-called dual observables arising from a vari- 
ation of the boundary conditions of the fermions in time- 
like direction have been introduced [42| and employed for 



L[A ] = — V exp 



dx A (x ,x) 



(1) 



1 Strictly speaking, we have to distinguish between the background 
temporal gauge field in Landau-DeWitt gauge and its expecta- 
tion value associated with the order parameter for confinement, 
see Refs. [33I. |39|| . We skip this subtlety here since it is of no im- 
portance for the present paper and refer to {Aq) as the position 
of the ground-state of the order-parameter potential. 
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Figure 1: Normalized order-parameter potential V(/3{Ao)) for SU(2) (lef t panel) and SU(3) (right panel) Yang-Mills theories 
for various temperatures as obtained from a first-principles RG study [331 ]. For SU(2) we show the potential for temperatures 
between T = 260 MeV and T = 285 MeV (from bottom to top). For 5(7(3) the Cartan subalgebra is two-dimensional and, in 
turn, the potential depends on two independent variables. Here, we show only a slice of the potential for various temperatures 
between T — 285 MeV and T = 310 MeV (from bottom to top) in the relevant direction of the Cartan subalgebra. The first- 
order phase transition in SU(3) Yang-Mills theory is indicated by a jump in the position of the minimum, see also the inlay of 
the figure. The phase transition temperature is T<j ~ 266 MeV for SU(2) and Td « 290 MeV for SU(3), respectively. 



a study of the relation of quark confinement and chiral 
symmetry breaking at finite temperature [3oT - [38l |43U50| . 
These dual observables relate the spectrum of the Dirac 
operator to the order parameter for confinement, namely 
the Polyakov loop. The introduction of these observables 
constitutes an important formal advance which allows us 
to gain a deeper insight into the underlying dynamics at 
the QCD phase boundary. However, they do not allow 
us to fully resolve the question regarding the relation of 
quark confinement and chiral symmetry breaking. 

In this work we aim to shed more light on the question 
under which circumstances the chiral and the deconfinc- 
ment transition lie close to each other 2 . To this end, we 
analyze the deformation of the RG fixed-point structure 
of chiral four-fermion interactions due to confining gauge 
dynamics in Sect. HU As detailed in Refs. (15l-Tl9l l5l| 
these four-fermion interactions can then be easily con- 
nected to the QCD Lagrangian at high momentum scales. 
In Sect. IIIII we then present our results for a partially 
bosonized formulation of the ansatz discussed in Sect. UT1 
In particular, a necessary condition for an exact mapping 
of both the fermionic and the partially bosonized theory 
is discussed. As one of our main results, we present a 
phase diagram spanned by the temperature and the pion 
decay constant f w . The latter is directly related to the 
chiral condensate. This phase diagram allows us to gain 
some insight into the interrelation of quark confinement 
and chiral symmetry breaking. In contrast to the phe- 
nomenologically more relevant phase diagram spanned 
by the temperature and quark chemical potential, the 



2 In full QCD with light but finite quark masses both transitions 
are crossovers. Since there is no unique way to define the criti- 
cal temperature associated with a crossover, a proof of an exact 
coincidence of the two transitions seems to be impossible in any 
case. 



(T, fir) phase diagram can be studied by various differ- 
ent approaches without suffering from problems, e. g., 
arising from a complex-valued Dirac operator. A study 
of the (T, f„) phase diagram may therefore be helpful to, 
e. g., benchmark continuum approaches with the aid of 
lattice simulations. Our concluding remarks and possi- 
ble future extensions of the present study are given in 
Sect. E| 



II. FERMIONIC FIXED-POINT STRUCTURE 
AND THE LOCKING MECHANISM 

We start our discussion of a dynamical locking mech- 
anism for the chiral phase transition with an analysis of 
the fixed-point structure in the matter sector of QCD. For 
our study we employ an RG equation for the quantum 
effective action, the Wetterich equation [52|. The effec- 
tive action T then depends on the RG scale k (infrared 
cutoff scale) which determines the RG 'time' t = ln(fc/A) 
with A being a UV cutoff scale. For reviews on and in- 
troductions to this functional RG approach we refer the 
reader to Refs. [53M6ll ] . 

For our more general discussion in this section, it suf- 
fices to consider the following ansatz for the effective ac- 
tion: 

T k [^,^,(A )} = J d 4 x{z^(ip+g l0 (A })^ 

+^ m? - (Hi^) 2 ] } , (3) 

For Xff, = this ansatz can be considered as the mi- 
croscopic matter part of the QCD action functional in 
Landau-DcWitt gauge. In the present work we restrict 
ourselves to Nf = 2 massless quark flavors and N c col- 
ors. The Tj's represent the Pauli matrices and couple 
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the spinors in flavor space. The coupling is consid- 
ered to be RG-scale dependent. Note that fermionic self- 
interactions are fluctuation-induced in full QCD, e. g. by 
two-gluon exchange, and are therefore not fundamental, 
see Refs. (l5l - [20| for a detailed discussion. 

Our ansatz for the effective action can be considered 
as the leading order in a systematic derivative expansion 
of the fermionic sector of QCD. The associated small pa- 
rameter of such an expansion is the anomalous dimension 
= — dt In Z^ of the fermion fields which has indeed 
been found to be small in earlier studies with vanishing 
background field (A ), sec Rcfs. [H, [H H H ■ In the 
following we therefore set the wave- function renormaliza- 
tion Z$ = 1. 

In this section we also drop a possible momentum de- 
pendence of the four-fermion coupling, A^(|p| <C k). This 
approximation does not permit a study of properties, 
such as the meson mass spectrum, in the chirally broken 
regime; for example, mesons manifest themselves as mo- 
mentum singularities in the four-fermion couplings. How- 
ever, the point-like limit can still be a reasonable approx- 
imation in the chirally symmetric regime above the chiral 
phase transition which allows us to gain some insight into 
the question how the theory approaches the regime with 
broken chiral symmetry in the ground state jl7H19{ . For 
our more quantitative analysis in Sect. lIIll we shall partly 
resolve the momentum dependence of the fermionic inter- 
actions in order to gain access to low-energy observables. 
At this point it is important to stress that in the point- 
like limit the RG flow of the four-fermion coupling, which 
signals the onset of chiral symmetry breaking, is com- 
pletely decoupled from the RG flow of fermionic n-point 
functions of higher order. For example, 8-fermion inter- 
actions do not contribute to the RG flow of the coupling 
A^, in this limit. 

The constant background field (Aq) in our ansatz §3§ 
is an element of the Cartan subalgebra and can be 
parametrized in terms of the corresponding generators 

f3g(A Q ) = 2* E T a ^ 

T a 6Cartan 

= 2n T a v (a) \<f>\, v 2 = l, (4) 

T a eCartan 

where the T a 's denote the generators of the underlying 
SU(N C ) gauge group in the fundamental representation 3 . 
It is convenient to introduce the eigenvalues v% of the 
hermitian matrix in Eq. (TJ|: 

n =spcc{(TV% | v 2 = 1} . (5) 

The presence of a finite background field (Aq) now yields 
a fermion propagator (r^ 2 )) -1 (inverse two-point func- 
tion) which is no longer proportional to the identity 1 in 



3 The dimension of the Cartan subalgebra is N c — 1. 




Figure 2: Representation of the 1PI Feynman diagram associ- 
ated with the A^-term on the right-hand side of the RG flow 
equation (01. Our functional RG study includes resumma- 
tions of this diagram to arbitrary order in A^. 

color space. However, it can be spanned by the Cartan 
subalgebra as follows: 

(r^K^I}])" 1 = ^ (r^ 2) K^I0l>]) 1 i« 

+ E (ri 2) [{^|0|}])" 1 ^,(6) 

T a e Cartan 

Here, the TjS's denote the generators in the fundamental 
representation. The expansion coefficients on the right- 
hand side of Eq. ([6]) can be obtained straightforwardly 
by using ti F T a T b = \8 ab and tr F T a = 0. Of course, it is 
expected from a physical point of view that the expan- 
sion (|6]) of the fermion propagator is convenient once a 
background field (Aq) is introduced into the theory since 
the latter distinguishes a direction in color space. For 
(A ) = wc have (r^)- 1 = 0. 

Using the ansatz ([3]) together with the parametriza- 
tion (||)]), we obtain the RG flow equation for the di- 
mensionless renormalized four-fermion coupling A^ in the 
point-like limit: 

= 81X4, = (2 + 2?7^)A V , 

-4( 2 +^)E^ F) (^O,^H)A 2 ,(7) 
n c i=i 

where 

A,/, = Z^j 2 k 2 \^ . (8) 

Note that A^, depends on the background field (Aq) and 
the dimensionless temperature r = T/k. The so-called 

threshold function l\ corresponds to a one-particle irre- 
ducible (1PI) Feynman diagram, see Fig. [21 and describes 
the decoupling of massive and thermal modes. Moreover, 
the regularization scheme dependence is encoded in these 
functions. The definition of the threshold function 
can be found in App. [A] 

Let us now discuss the fixed-point structure of the cou- 
pling A^. Apart from a GauBian fixed point we have a 
second non-trivial fixed point. The value of this fixed 
point depends on the dimensionless temperature r and 
the dimensionless coordinates of the background 
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field (Aq), see Fig. [3l At vanishing temperature (and 
background field (Aq)) we find 



A* 



(2JV c + l)lf (0,0,0) 



(2A C 



1] 



(9) 



for the non-Gaufiian fixed point. For illustration we have 

evaluated the threshold function l\ in the second line 
for the regulator function (|A2[) . Here, rjl denotes the 
value of the fermionic anomalous dimension at the fixed 
point. Note that the rescaled fixed-point coupling iV c A3. 
approaches a constant value in the limit N c — > oo. 

The fixed-point value is not a universal quantity as 
its dependence on the threshold function indicates. How- 
ever, the statement about the existence of the fixed point 
is universal. Choosing an initial value A^[ v < A^, at the 
initial UV scale A we find that the theory becomes non- 
interacting in the infrared regime (A^ — > for k — > 0), 
see Fig. [3] For A^ v > A^ we find that the four-fermion 
coupling A^, increases rapidly and diverges eventually at 
a finite scale fe cr . This behavior indicates the onset of 
chiral symmetry breaking associated with the formation 
of a quark condensate. Hence chiral symmetry break- 
ing in the IR only occurs if we choose A^ v > A3,- Of 
course, the divergence of the four-fermion coupling at a 
finite scale fc cr is an artifact of our point-like approxima- 
tion. It can be resolved by taking into account (some 
of) the momentum dependence of the coupling A^, see 
Sect. IIII1 This will then allow us to gain access to QCD 
low-energy observables. In any case, the scale fc cr at 
which I/A^,(fc cr ) = sets the scale for a given IR ob- 
servable O: 



O ~ k d ° , 



(10) 



where do is the canonical mass dimension of the observ- 
able O. At vanishing temperature the scale fc cr can be 
computed analytically. We find 

/ ^uv _ y \ I 

fc cr = Ae(Ag v -A;)( v '^ v t \ +o( v ;). (a) 

Thus, the critical value fc cr scales with the distance of 
the initial value A^ v from the fixed-point value A^,. For 
increasing A^ v the scale fc cr increases and, in turn, the 
values of low-energy observables, such as the pion de- 
cay constant f v and the chiral phase transition temper- 
ature T x , increase. 

From now on, we assume that we fix A^ v > A^ at 
T = 0. The value of A^J V then determines the scale 
fccr = fccr(A^ v ) which is related to the values of the low- 
energy values. For a study of the effects of a finite tem- 
perature and a finite background field (A ) we then leave 
our choice for A^ v unchanged. This ensures comparabil- 
ity of the results at zero and finite temperature for a given 
theory defined by the choice for A^J V at zero temperature. 




(A ) > 



Figure 3: Sketch of the /3 function of the four-fermion in- 
teraction for three different cases: vanishing temperature 
(black/solid line), a given finite value of the temperature T 
and (Ao) = (red/dashed line), the same temperature T but 
(Ao) > (blue/dashed-dotted line). The arrows indicate the 
direction of the RG flow towards the infrared. 



Next, we turn to a discussion of the fixed-point struc- 
ture at finite temperature but vanishing background 
field (Ao) . We still have a Gaufiian fixed point. Moreover, 
we find a pseudo fixed-point A3,(t) for arbitrary values 
of t at which the right-hand side of the flow equation is 
zero: 



K(r) 



(2A c + l)4 F) (r,0,0) 



(12) 



For high temperatures T k wc find A3, ~ (T/k) 3 . Let 
us now assume that we have chosen A^ v > (r = 0) . 
Since the the value of the (pseudo) fixed-point increases 
with increasing T/k, the rapid increase of the four- 
fermion coupling towards the IR (k — > 0) is effectively 
slowed down and may even change its direction on the 
plane defined by the two-dimensional /3\ , function 4 , see 
also Fig. |3l This behavior of the pseudo fixed-point (t) 
already suggests that for a fixed initial value A^ v a crit- 
ical temperature T x exists above which the four-fermion 
coupling docs not diverge but approaches zero in the IR. 
Such a behavior is indeed expected for high temperatures 
since the quarks become effectively stiff degrees of free- 
dom due to their thermal mass ~ T and chiral symmetry 
is restored. 

Let us now turn to a discussion of the fixed-point struc- 
ture for finite T and (Aq). In the present approach we 
consider the value of the background field (Aq) as an 
external input which is determined by the ground state 
of the corresponding order-parameter potential, see also 
Fig. [TJ As discussed above, the position (Aq) of the 
ground-state is directly related to our order parameter for 
confinement, trF^[(Ao)]. For temperatures much larger 
than the deconfincment phase-transition temperature T<j 
we have (Aq) = 0, i. c. trp-k^o)] = 1- On the other 



At finite temperature the f}\ 
namely r = T/k and A^,. 



function depends on two variables, 
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hand, the position (Aq) of the ground state in the con- 
fined phase of pure SU (N c ) Yang-Mills theory is uniquely 
determined up to center transformations by [lU H^] 

tr F (L[(A )]) ra = (13) 

with (n mod N c ) = 1, , . , , N c — 1. These conditions de- 
termine the N c —1 coordinates {0^} of (Aq), see Eq. (jH). 
Moreover, we have 

-l|tr F (i[(A )]r|<4- (W) 



for n € N. In QCD with dynamical quarks the ground- 
state value (Ao) is then shifted and yields a small but 
finite order parameter in the confined phase. 

For temperatures much larger than the deconfinement 
phase-transition temperature Td, the fixed-point struc- 
ture remains unchanged since (.Ao) tends to zero for 
T > T d . For finite (Aq) the pseudo fixed-point A^ de- 
pends on (Aq) and r = T/k. The position of the pseudo- 
fixed point AJ can again be given in closed form: 



A?, (r, (/!„)) 



^ 2 4)| ,sF,M 'H 



1 / 1 \ 00 

^2 2 + E(-^)" [M£[<4>)]) B + tr P (£t[(A )])"] (l + J) e"? 



A,* (0,0) 6tt2 V AT, 



. (15) 



where we have dropped terms depending on ryj, on the 
right-hand side. To obtain the second line we have em- 
ployed the regulator function (|A2[) . However, the general 
form of the asymptotic series (|15p holds for any regula- 
tor function as can be shown by means of Poisson re- 
summation techniques 5 . Note that the series (|15[) can be 
considered as an expansion for small r = T/k. 

Using Eq. (fT3"|) we observe that all finite-temperature 
corrections to the fixed-point value vanish identically 
in the confined phase for N c — > oo, provided that the 
ground-state value (Ao) is identical in SU(N C ) Yang- 
Mills theory and QCD with dynamical fermions. Of 
course, the latter assumption is not exactly fulfilled but 
for physical quark masses it is reasonable to assume 

tr F L[(A )] « 1 (16) 

at low and intermediate temperatures. Thus, we have 
found that 

a;(o,o) ee x;(t,(a )) (17) 

in the limit N c — > oo, independent of the temperature T 
for T < Td- On the other hand, we have 

A;(t,(A »^A;(t,0) for {Ao)->0. (18) 

With the same reasoning we find 

/3 a ,(0,0)ee/3^(t,(A )) (19) 

for T < T d and N c — > oo. This means that for T < T d 
the question of whether chiral symmetry is spontaneously 



5 The sum over the r-dependent terms in the second line of 
Eq. I I15I I is also closely related to the geometric series. 



broken or not is in fact independent of the temperature, 
but depends only on the choice of the initial value A^ v 
relative to its fixed-point value AJ at T = and (Aq) = 0. 
We emphasize that Eqs. (fT7)) - (|T5]) arc rcgularization- 
schemc independent statements. For any admissible so- 
called 3d regulator function, e. g. Eq. (|A2|> . it is also 
possible to show analytically that 

a;(o,o)<a;(t,(A)))<a;(t,o) (20) 

for arbitrary temperatures r = T/k for N c — 2, 3; see 
also discussion of Eq. (|A10[) in App. [5] In any case, 
X^(t,(Ao)) interpolates continuously for a given finite 
value of r between AJ(0,0) and X1(t,0). 

Provided that we choose an initial value A^[ v > 
A^(0, 0), it follows immediately from Eq. (fT9j) that 

T x > T d (21) 

for N c —> oo, see also Ref. 0]. This means that the 
chiral phase transition is locked in due to the confining 
dynamics in the gauge sector. Loosely speaking, thermal 
fluctuations of the quark fields, which tend to restore 
chiral symmetry, are suppressed since they are directly 
linked to the deconfinement order parameter. Thus, we 
have found that the restoration of chiral symmetry is 
intimately connected to the confining dynamics in the 
gauge sector. 

In the present analysis the actual chiral phase tran- 
sition temperature depends on two parameters, namely 
the value of the background field (Ao) and the initial 
condition A^J V . However, Eq. (|2Tj) is a parameter- free 
statement. It simply follows from an analysis of the ef- 
fect of gauge dynamics on the fixed-point structure in the 
fcrmionic sector. In particular, we have only made use of 
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general properties of the deconfinement order parameter 
and the fact that A,^ v > A^(0, 0) is a necessary condition 
for chiral symmetry breaking at T = and (Aq) = 0. Of 
course, the initial condition A^ v is not a free parameter 
in QCD but originally generated by quark-gluon interac- 
tions at high (momentum) scales. In a given rcgulariza- 
tion scheme the value of A^ v can therefore in principle 
be related to the value of the strong coupling a s at, e. g., 
the r mass scale (l6l - fl9j . We would like to point out that 
neither the value of a s at some scale nor the value of A^J V 
on a given RG trajectory is a physical observable. How- 
ever, their values can be related to physical low-energy 
observables. Recall that the value of A^ v determines the 
critical scale fc cr which sets the scale for IR observables, 
see Eq. (|TT]l . 

Our findings in the limit N c — > oo even allow us to es- 
timate a window in parameter space in which the chiral 
phase transition and the deconfinement phase transition 
lie close to each other. In our discussion of chiral symme- 
try breaking for vanishing background field (Aq) we have 
argued that T x ~ fc cr , where the scale k cr is eventually 
determined by our choice for A^ v . If we take into account 
the background field (Aq), then the chiral phase transi- 
tion temperature T x is locked in and we necessarily have 
T x > T d in the large- N c limit, see Eq. ([2l]), Thus, the 
chiral phase transition temperature for all theories which 
would allow for T x < Td for (Aq) = is shifted such that 
T x ~ Td. The upper end of the locking window can there- 
fore be estimated by the smallest value for A^ v for which 
T x for vanishing (Aq) is still larger than Td. Whereas 
Al and A^ v are scheme-dependent quantities, the mere 
existence of such a window in parameter space is a uni- 
versal statement. As we have argued above, the value 
A^ v can be related to the values of physical low-energy 
observables. Therefore the existence of such a window 
for the initial condition A^ v suggests the existence of a 
corresponding window for the values of low-energy ob- 
servables. We shall come back to this below when we 
discuss the phase diagram in the plane spanned by the 
temperature and the pion decay constant. 

Let us now discuss the relation of quark confinement 
and chiral symmetry breaking for finite N c . In this case, 
terms with 

n mod N c = (22) 

contribute to the right-hand side of Eq. (| 15[) and to the 
RG flow of A^,. We then find that the inequality (j20|) 
holds only for r = T/k <C 1 but not for arbitrary values 
of t. However, this does not necessarily imply that we 
do not have a finite range of values for the initial condi- 
tion A^J V anymore in which the chiral and the deconfine- 
ment phase transition are tightly linked. It only implies 
that the lower end of the window for A^ v is shifted to 
larger values compared to the large- N c limit. 

To illustrate our analytic findings we have studied nu- 
merically the RG flow of the four-fermion coupling A^ for 
finite N c . In our setup, the phase transition temperature 
is defined to be the smallest temperature for which A^, 
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Figure 4: Phase diagram for two massless quark flavors and 
N c — 2 colors (blue/dashed line) as well as for N c = 3 colors 
(black/solid line) in the plane spanned by the temperature 
and the rescaled coupling A^ v /A^,. The lines depict our re- 
sults for the ratio T x /T^ of the chiral and the deconfinement 
phase transition temperature as a function of A^ v /A^,. Recall 
that there is no splitting of the phase boundary (i. e. T x ~ Td) 
for small A^ v in the limit N c — ¥ oo. 



remains finite in the infrared limit k — > 0. Strictly speak- 
ing, this only yields an upper bound for the phase transi- 
tion temperature since it is only sensitive to an emergence 
of a condensate on intermediate momentum scales but in- 
sensitive to a fate of the condensate in the deep IR due 
to fluctuations of the Goldstone modes In Fig. 0] we 
present the phase diagram for two massless quark flavors 
and N c = 2 as well as N c = 3 in the plane spanned by the 
temperature and the UV coupling A^ v . To estimate the 
phase boundary we have employed the results for (Aq) for 
the corresponding SU(N C ) Yang-Mills theory as obtained 
in Refs. [33|, H^], see also Fig. [TJ The associated decon- 
finement phase transition temperature is Td ~ 266 MeV 
for pure gauge ST/ (2) and Td ~ 290 MeV for pure gauge 
SU(3), respectively. For the UV cutoff we have chosen 
A = 1 GcV. In accordance with our analytic findings 
we observe that there is only a chirally symmetric phase 
for A^J V < Al. Taking into account corrections beyond 
the large- N c approximation, we find that the chiral phase 
transition temperature increases continuously with A^ v , 

T x^(A" V -A;,) 1/2 , starting from T x = for Aj^ v = AJ. 
In this regime we have Td > T x . Increasing AyJ v further 
we find a window for the values of Aj^ v in which we have 
T x ~ Td. In this regime the chiral phase transition is 
locked in due to the confining dynamics in the gauge sec- 
tor. Note that T x is a strictly monotonously increasing 
function in A^ v for (Aq) =0. In agreement with our 
analytic results we observe that the size of this locking 
window increases with increasing N c . For A^ v 3> XI we 
then observe that the phase transition temperatures Td 
and T x differ again. However, we now have T x > Td. 

At this point a few comments are in order. First, it is 
clear that the potential of the confinement order param- 
eter in full QCD (and hence the position of the ground- 
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state (Aq)) also receives contributions from Feynman dia- 
grams with at least one internal fermion line. These con- 
tributions tend to lower the deconfinement phase tran- 
sition temperature 6 . We stress that we only use the 
Yang-Mills approximation for (Aq) here to explore the 
impact of the discussed locking mechanism for the chiral 
phase transition on the finite-temperature phase struc- 
ture. In any case, our analytic findings are not (strongly) 
affected by this approximation since they rely on very 
general properties of the confinement order parameter. 
Therefore we still expect that a window in parameter 
space exists in which the chiral and the deconfinement 
phase transition lie close to each other. However, our 
estimate of the phase diagram in the (T, A^ v )-plane and 
the size of the locking window will change quantitatively 
when we take into account the corrections to (Aq) due 
to quark fluctuations. In particular, we expect that be- 
low the locking window the difference T x and Td be- 
comes smaller since this regime is associated with quarks 
with a small dynamical mass. Second, our ansatz ([3]) for 
the effective action is not complete with respect to Fierz 
transformations, see e. g. Refs. H3] ; for example, 

we have dropped the so-called vector-channel interaction 
~ (Tpj^ip) 2 . Such interactions would also contribute to 
the RG flow of the four-fermion interaction A,/, . At finite 
temperature the minimal basis of point-like four-fermion 
interactions is larger than at vanishing temperature since 
the Poincare invariance is broken by the heat bath. If we 
allow for a finite (Aq), the minimal set of point-like four- 
fermion interactions is even larger than in the case of 
vanishing background field (Aq). This is due to the fact 
that a finite background field (Aq) distinguishes a di- 
rection in color space. For example, our expansion of 
the fermion propagator suggests that a finite background 
field (Aq) gives rise to additional point-like interactions of 
the type - (^T( 3 V) 2 and - (V?T( 8 V) 2 for N c = 3. How- 
ever, the additional diagrams arc of the same topology as 
the one shown in Fig.[2J We therefore expect that the in- 
clusion of additional four-fermion interactions associated 
with a Fierz-complete basis is particularly important for 
a quantitative computation of the chiral phase transition 
temperature beyond the large- N c limit. Such an analysis 
is beyond the scope of this work and deferred to future 
studies. 

Finally we would like to briefly comment on 
PNJL/PQM-type model studies which are closely related 
to the present study. In contrast to the present work, 
the PNJL/PQM-type model studies use the assump- 
tions tr F L[(A )} = (tr F L[A }) and N?tr F (L[(A )]) n = 
N c (trpL[Ao]) n , see e. g. Ref. @. If we used these as- 
sumptions, we would not have observed the regime with 
T x < Td for small values of (A^ v — Al) and finite N c where 
the two phase transitions are decoupled. As we have dis- 



Strictly speaking we only have only a deconfinement crossover in 
the presence of light quarks. 



cussed above, the decoupling of the phase transitions for 
small values of (A^ v — AJ) is only absent in the limit 
A c — > oo in our study. Therefore we conclude that these 
approximations used in PNJL/PQM-type model studies 
correspond to a large- A c approximation in the coupling 
of the matter and the gauge sector. This type of large- N c 
approximation should not be confused with the large- N c 
approximations used in the matter sector of these mod- 
els, such as neglecting pion fluctuations. In any case, 
our analysis shows that the large- A c approximation as- 
sociated with the approximation trpi[(Ao)] = (trpi[Ao]) 
clearly affects the dynamics near the finite-temperature 
phase boundary and may therefore also affect the predic- 
tions for the (T, n) phase diagram from PNJL/PQM-type 
models. 



III. THE LOCKING MECHANISM AND THE 
(T, f n ) PHASE DIAGRAM 

In the previous section we have discussed the interrela- 
tion of quark confinement and chiral symmetry breaking 
in a purely fcrmionic language. We have found that there 
exists a regime in the phase diagram spanned by the tem- 
perature and the coupling A,^ v in which the chiral and 
the deconfinement phase transition are tightly linked. In 
this section, we would like to map this phase diagram 
onto a phase diagram spanned by the temperature and 
an IR observable of QCD, e. g. the chiral condensate 
K^ip)] 1 / 3 or the pion decay constant / w . 

In QCD we have only one input parameter, e. g. 
the strong coupling a s at a (high) momentum scale or, 
cquivalently, Aqcd- For vanishing current quark masses 
Aqcd determines all physical observables, such as the 
pion decay constant f v ~ Aqcd as well as the decon- 
finement and the chiral phase transition temperature, 
Td ~ Aqcd and T x ~ Aqcd- Thus, real QCD in a 
phase diagram spanned by the temperature T/Aqcd and 
/tt/Aqcd is a single point. In the present paper the scale 
Aqcd is fixed by the input for the background field (^4o)- 
In contradistinction to real QCD, however, an additional 
parameter is present in our study, namely A^[ v . In the 
following we shall consider this parameter as an asset 
which allows us to deform QCD. In our model, two-flavor 
QCD as defined by f n « 90MeV then corresponds to a 
specific choice for A5J on a given RG trajectory. Thus, 
our deformed model effectively depends on two param- 
eters, namely A^J V and Aqcd- In particular, the chiral 
phase transition temperature depends on two parameters 
in our study, namely A^J V and Aqcd , whereas the decon- 
finement phase transition temperature depends only on 
Aqcd ~ Ti. In the remainder of the paper we exploit 
the dependence of our model on A^J V in more detail to 
gain some insight into the relation of chiral and confining 
dynamics close to the finite-temperature phase bound- 
ary. Eventually, this leads to a prediction for a (T, f^) 
phase diagram. We would like to add that deformations 
of QCD-like theories with an additional relevant parame- 
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ter, such as a four-fcrmion coupling, indeed play a promi- 
nent role in beyond-standard model applications, see e. g. 
Ref. |e3. 

For the computation of the (T, f n ) phase diagram we 
employ a (partially) bosonized version of our ansatz (|3|). 
Whereas the purely fermionic description is particularly 
convenient for analytic studies, the partially bosonized 
ansatz allows us to resolve momentum dependences of 
fermionic self-interactions in a simple manner and there- 
fore permits a study of the formation of the chiral con- 
densate and the mass spectrum of mesons. To obtain 
the partially bosonized formulation of our ansatz ([3]) we 
perform a Hubbard-Stratonovich transformation of the 
underlying path- integral [66l . Wi\ . This introduces auxil- 
iary fields <i T = (<t, 7?) into the theory which mediate the 
interaction between the quarks. Here, we assume that 
the bosons are composites of fcrmions and do not carry 
an internal charge, e. g. color or flavor: a ~ (ipip) and 
7? ~ (ipT^sip). The components of $ are labeled accord- 
ing to the role that the corresponding fields are playing 
in the spontaneously broken regime. The effective action 
of the partially bosonized theory then reads 7 

T k $,il>,$,{Ao)] = J d 4 x{z,,pi>(i@ + g la {A })^ 

1 - 2 - - 

+ -Z$ (c^$) + ihip(a+ir ■ 7?7 5 )V 

+ ^™ 2|>2 + ^** 4 }' (23) 

with a Yukawa coupling and the boundary condi- 

tions 



lim Z& 

fe^A 

lim Z^h 
lim A* 

fc^A 



0. 

1. 

0. 



(24) 
(25) 
(26) 



These boundary conditions together with the identity 



(27) 



allow us to map the ansatz (|23p onto the model (O at 
the initial UV scale A. In particular, we are now able to 
relate the initial values A^ v to physically meaningful IR 
observables, e. g. the pion decay constant and the chiral 
condensate. 



The term ~ $ 4 in the effective action (|23|) corresponds 
to an 8-fermion interaction term in the purely fermionic 
description, <£> ~ "tpip. Due to the boundary condition (|26|) 
this term is generated dynamically and not ad j uste d by 
hand in our RG approach, see e. g. EI 0, Hi M if. 
Thus, the value of the corresponding coupling at the 
initial RG scale A does not represent an additional pa- 
rameter of the theory. In any case, this coupling only 
plays a prominent role in the deep IR regime where it ac- 
counts for the mass difference between the c-meson and 
the pions. 

A word of caution concerning the mapping of the 
partially bosonized and the purely fermionic description 
needs to be added here. The identity ([27|) suggests that 
we indeed have only one input parameter in the present 
study, namely the value of the ratio h^y/m^y at the UV 
scale A. In d = 4 space-time dimensions, however, the 
Yukawa coupling h is marginal. This suggests that the 
partially bosonized theory in d = 4 depends on two in- 
put parameters in contrast to d = 3, see Refs. (o3l. lo8l]. 
Nonetheless the critical scale k CI only depends on the 
ratio ^uv/™uv m leading order in an expansion in pow- 
ers of 1/N C and receives only small corrections from the 
next-to-leading order, see below. On the other hand, the 
ratio of IR observables, such as the ratio of the a mass 
and the constituent quark mass, depends on both param- 
eters [iH. In the following we choose the ratio ^uv/^uv 
and the Yukawa coupling h\jy as the independent input 
parameters at the UV scale A. 

Let us now discuss the RG flow equations of the di- 
mensionless renormalized couplings m 2 = m 2 /(Z$/c 2 ), 
A$ = \^,jZ\ and the Yukawa coupling h = h/ {z]J 2 Z^). 
Using the ansatz (|23|) together with the parametriza- 
tion (JB|) of the fermion propagator we find 



dtm 2 



3 

(i]<j> - 2)m 2 - — -rZi(r, m 2 )A$ 
2ir z 



N, 



<9 t A$ = 



d t h 2 



i=i 
3 



2?7$A$ + -^h{r, m 2 )A| 
(=1 

(77$ + 2rj^)h 2 
2 1 N < 



(28) 



(29) 



^£Ci B W>*M,m> 4 ,(30) 



1=1 



7 At finite temperature the Poincare invariance of the theory is 
broken. Therefore the wave-function renormalizations longitudi- 
nal and transversal to the heat-bath obey a different RG running. 
We neglect this difference in the present work. This is justified 
since it has indeed been found in Ref. [63l | that the difference is 
small at low temperatures and only yields mild corrections to, 
e. g., the thermal mass of the bosonic degrees of freedom for 
intermediate temperatures T >T X . 



where 77$ = — i9 t lnZ$ and rj^, = — 9 t lnZ^. The thresh- 
old functions are defined in App.|XJ Note that the quarks 
are massless in the symmetric regime, m q = 0. 

At this point it is instructive to have a closer look 
at the mapping between the partially bosonized and the 
purely fermionic description. To this end, we consider 
the RG flow of the ratio h 2 /m 2 which can be obtained 
straightforwardly from the flow equations (f2"8")l and (|30|) . 
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Wc obtain 



h 2 \ ,_ _ , f h 2 \ 3 . , 2 , . ( h 2 



d t[ — ) = (2 + 2 %) ( )+^h(T,m')\* [ — 



rrr 

x 2 



3FEff( f X.*™ > )(S)-( 31 ) 



7T 2 N. 



1=1 



Using Eqs. (US]) and {57) as well as 



(FB) 
1,1 



(m>l) 1 ( F ) 



2 l l 



(T,m>#|), (32) 



we recover the RG flow equation ([7]) of the four-fermion 
coupling A,/,. Thus, the partially bosonized and the 
purely fermionic description are indeed identical at the 
UV scale A. Note that the prefactor of the term ~ h 2 /m 2 
would turn out to be incorrect if we did not include the 
RG running of the Yukawa coupling. In other words, a 
standard local potential approximation does not incorpo- 
rate all terms associated with a systematic expansion of 
the flow equations in powers ofl/N c . In Ref. J68| this ob- 
servation is discussed in the context of the Gross-Neveu 
model. 

From Eq. wc also deduce that the partially 

bosonized description allows us to go conveniently be- 
yond the point-like approximation employed in the purely 
fermionic description discussed in Sect. HU To be more 
specific, we observe that the momentum dependence of 
the four-fermion vertex is effectively parametrized by the 
RG flow of the four-boson coupling A* , the Yukawa cou- 
pling h and the mass parameter m 2 . We stress that the 
purely fermionic point-like description and the partially 
bosonized description are no longer identical for scales 
k < A. From the flow equation of the ratio h 2 /m 2 , 
however, it follows that the differences are quantitatively 
small in the symmetric regime since the renormalized 
mass of the mesons is large. 

In the regime with broken chiral symmetry in the 
ground state the mass parameter m 2 assumes negative 
values. This behavior signals the existence of a finite vac- 
uum expectation value of the field In this regime, it is 
therefore convenient to study the RG flow of the four- 
boson coupling A$ and the vacuum expectation value 
($) = $ . The RG flow of the latter can be obtained 
from the following condition: 



d_ 



0. 



(33) 



$0 



The flow equations for $ an d A$ are then given by 



Av 



^l^^ml^l)—, (34) 
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Figure 5: Phase diagram for two massless quark flavors and 
iV c = 3 in the plane spanned by the rescaled temperature 
T x /Td and the value of the pion decay constant f n at T = 0. 
The line depicts our result for the ratio T x /Ta of the chiral 
and the deconfinement phase transition temperature. 



9 3 
<9 t A$ = 2?7$A$ + — -h(T, ml)X% + -~-^h{r, m 2 )X% 



1=1 



where 



= A$$q, ml = and m 2 = /i 2 $q is the con- 

— 1/2 - 

stitucnt quark mass. Note that we identify <I>o with 

the pion decay constant For simplicity, we neglect 
the running of the Yukawa coupling in the regime with a 
broken chiral symmetry of the ground-state. In fact, 1PI 
diagrams with at least one internal quark line are para- 
metrically suppressed in this regime since the quarks ac- 
quire a (large) mass. Therefore the RG flow of the theory 
in the spontaneously broken regime is mainly governed 
by 1PI diagrams with internal pion lines only. However, 
the latter class of diagrams does not directly contribute 
to the RG flow of the Yukawa coupling 8 . 

In the present study we also neglect the running of the 
bosonic and fermionic wave- function renormalizations. 
Thus, we consider Z$ = I and = 1. In studies with 
vanishing gluonic background field [HI, [T^, [62j it has in- 
deed been found that the associated anomalous dimen- 
sions 77$ and rjtp are small over a wide range of scales. 
Nonetheless using Z$ = 1 means that we violate the 
initial condition m 2 3> 1 for the renormalized mass at 
the UV scale A. The latter condition is equivalent to 
the condition (f2"4|) . Therefore our study of the partially 
bosonized theory necessarily relies on two parameters in- 
stead of one, e. g. the initial values for h 2 and h 2 /m 2 , 
which wc expect to be anyway the case, see discussion 



!=1 



8 In the symmetric regime we have only taken into account the 
contributions to the running of the Yukawa coupling h which are 
required to map the flow of the ratio h? /m 2 onto the RG flow of 
the four-fermion coupling A^, at the UV scale A. 
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above. Employing Z$ = 1 then seems to be a reasonable 
approximation for a first explorative study of the locking 
mechanism discussed in Sect. HU In contrast to the flow 
of Z$, the RG flow of Z$ is solely driven by 1PI dia- 
grams with at least one internal boson and one internal 
fermion line. Since the mesons are heavy in the symmet- 
ric regime and the fermions are heavy in the regime with 
broken chiral symmetry in the ground state, it is justified 
to consider = 1 for our purposes. For a more quan- 
titative study including a determination of critical expo- 
nents, however, the running of Z^, Z$ and the Yukawa 
coupling needs to be taken into account in both regimes 
since very close to the phase transition both the mesons 
and the quarks are approximately massless. The corre- 
sponding RG flow equations can be derived along the 
lines of, e. g., Refs. [ll M, HI • 

Using the partially bosonized theory we now estimate 
the size of the locking window for the chiral phase tran- 
sition in terms of the pion decay constant The lat- 
ter is directly related to the chiral condensate | ("ipip) | 1//3 
via the Gell-Mann-Oakes-Renner relation. In Fig. [5] we 
present our result for the phase boundary for two mass- 
less quark flavors and N c = 3 in the plane spanned by 
the temperature and the value of the pion decay constant 
at T = 0. This phase diagram corresponds to the phase 
diagram discussed in Sect. [TTJ To obtain Fig. [5] we have 
set A = 1 GcV and used /itjv = 3 as initial condition 
at k = A. Different initial values for A^ v = /iuv/ m uv 
then translate into different values for f^. In addition, 
we have employed the results for (Ao) as obtained from 
a RG study of SU(3) Yang-Mills theory [H H| which 
is in very good agreement with lattice QCD results. We 
neglect the back-reaction of the matter sector on (Ao) 
and use only the Yang-Mills approximation for (Ao) here 
to explore the impact of the discussed locking mechanism 
on the finite-temperature phase structure. 

In Fig. [5] we notice that the chiral phase transition 
temperature tends to zero for — > and that T x is 
a monotonic function for small values of /„■, T x ~ f n . 
For < U < 100 MeV and /„ > 150 MeV we ob- 
serve a splitting of the phase boundary. Interestingly, 
we have T x > T d for f n > 150 MeV and T x < T d for 
fir ^ 100 MeV. As a non-trivial result we find that the 
chiral phase transition is locked to the deconfinement 
phase transition for 100 MeV < f n < 150 MeV. 9 This 
observation confirms our results from the study of the 
fixed-point structure in the purely fermionic description. 
Moreover, we find that the physically relevant range of 
values of the pion decay constant lies right below the 
locking window. The near coincidence of the two phase 
transition temperatures for ~ 90 MeV has also been 



Using h\jy = 3.5 as initial condition for the Yukawa coupling, 
we find that the value of the lower end of this window for /„■ is 
lowered by approximately 10%. However, the size of this window 
is roughly independent of the considered initial conditions for the 
Yukawa coupling. 



found in an RG study including the back-reaction of the 
quark fluctuations on the order-parameter potential j37j . 
Recall that for f n < 150 MeV the chiral phase transi- 
tion temperature is shifted to larger values compared to 
studies with (Aq) = since the thermal excitation of 
quarks is suppressed by the confining dynamics in the 
gauge sector, see Eq. (jT5j) . Loosely speaking, we observe 
a competition between the confining dynamics and the 
thermal excitation of the system. While an increase of 
the temperature tends to restore chiral symmetry, the 
confining dynamics in the gauge sector favors a ground 
state with broken chiral symmetry. 

F° r Jtt 150 MeV, which corresponds to constituent 
quark masses m q > 450 MeV, we find a second-order chi- 
ral phase transition for vanishing current quark masses. 
Moreover, our result for T x as a function of f n does 
not depend on the presence of the gluonic background 
field (Ao) any more. In this regime, we expect that our re- 
sults are least affected by our approximation of neglecting 
the back-reaction of the matter sector on the deconfine- 
ment order parameter. Here, the dynamics in the gauge 
sector and the matter sector are expected to be effec- 
tively decoupled. Nonetheless the difference between T x 
and Td might be smaller in a full QCD study. 

For f-n ^ 150 MeV a window may open up in which 
the chiral phase transition is of first order. This is due 
to the significant shift of the chiral phase transition tem- 
perature to larger values in this regime compared to the 
case without a gluonic background field. Provided the 
deconfinement order parameter rises quickly above the 
(deconfinement) phase transition [H, [H, [H, |69[ , strong 
thermal excitations may induce a chiral transition of first 
order. In our present study we indeed find a first-order 
chiral phase-transition for 100 MeV < f„ < 150 MeV. 
However, we rush to add that this observation is a short- 
coming of our approximations since we have employed 
the results for (A ) as obtained from a study of SU(3) 
Yang-Mills theory but neglected the back-reaction of the 
matter sector on (Aq). Thus, the first-order phase transi- 
tion in the gauge sector induces a first-order chiral phase 
transition in this regime. Nevertheless a first-order chi- 
ral phase transition may occur in this regime even if we 
include the back-reaction of the matter sector on the 
confinement order parameter, provided the latter rises 
rapidly for T > Td- For example, this might be the case 
for physical pion masses. In this respect the present anal- 
ysis provides a simple mechanism for a first-order chiral 
phase transition [70]. For f n < 100 MeV we find again a 
second-order chiral phase transition for vanishing current 
quark masses. 

We would like to remind the reader of the fact that we 
have only one parameter in QCD with massless current 
quarks and that four-fermion interactions are induced by 
fluctuations, e. g. two-gluon exchange; see Refs. h5l- 
r20l ] for a detailed discussion. In the present paper we 
have considered the initial value A^ v of the four-fermion 
coupling as an additional parameter which allowed us to 
deform QCD and study some aspects of the relation of 
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quark confinement and chiral symmetry breaking. In any 
case, we would like to stress that the relation of the fixed- 
point structure in the matter sector to the deconfinement 
order parameter is a universal statement which is a pri- 
ori independent of the details of the actual scale fixing 
procedure. This interplay of the matter and the gauge 
sector suggests the existence of a window for the values 
of the physical observables in which the chiral and the 
deconfinement phase transition lie close to each other. 



IV. CONCLUSIONS AND OUTLOOK 

In the present paper we have analyzed the interplay 
of the deconfinement and the chiral phase transition. To 
this end, we have studied the fixed-point structure of 
four-fermion interactions in the presence of a finite tem- 
poral gluonic background field. The latter can be di- 
rectly related to an order parameter for the deconfine- 
ment phase transition, namely trpL[(Ao)] > (trpi[v4o])- 
In a purely fermionic description the onset of chiral sym- 
metry breaking is indicated by rapidly increasing four- 
fermion interactions. Restoration of chiral symmetry is 
then displayed by four-fermion couplings which approach 
a Gaufiian fixed point in the IR. Thus, the question of chi- 
ral symmetry breaking in the IR can ultimately be linked 
to the fixed-point structure of the four-fermion couplings. 

We have indeed found that the fixed-point structure 
of four-fermion interactions is directly related to the or- 
der parameter for confinement. In particular, we have 
analyzed the scalar-pscudoscalar interaction channel and 
have shown that its (3 function docs not depend on the 
temperature in the limit N c — > oo, provided we are 
in the confined phase, i. c. ti'F£[(^4o)] ~> 0. Loosely 
speaking, thermal excitations of the quarks are sup- 
pressed in the confined phase. Therefore the chiral 
phase transition is locked in and shifted to higher tem- 
peratures. Our findings confirm the results of a mean- 
field study by Meisingcr and Ogilvie Q in which it has 
been pointed out that the chiral order-parameter po- 
tential is independent of the temperature in the con- 
fined phase by employing the assumption trp£[(Ao)] = 
(trpL[Ao]). Our non-perturbative analysis of the fixed- 
point structure of the matter sector establishes this ob- 
servation. However, we have not made use of the as- 
sumption tv-pL^Ao)] = (trpLLAo])- In fact, our study 
shows that the latter assumption is only justified in the 
limit N c — > oo. Moreover, we have studied how correc- 
tions to the large- N c approximation affect the dynamics 
at the finite-temperature phase boundary. The locking 
mechanism discussed in Sect. [TT] eventually allows us to 
determine a window for the values of the pion decay con- 
stant in which the chiral and the deconfinement phase 
transition lie close to each other. The size of this window 
depends on the number of colors N c and is maximal in 
the limit N c — > oo. 

For N c — 3 and two masslcss quark flavors we have 
computed the phase diagram in the plane spanned by the 



temperature and the pion decay constant. This phase di- 
agram can be divided into three different regimes. For 
small values of / T we find a regime with Xd > T x . For 
100 MeV < f„ < 150 MeV the deconfinement and the 
chiral phase transition lie close to each other. Finally, 
there is a third regime for > 150 MeV which is char- 
acterized by T x > Td- Here, the chiral condensate is 
large and therefore the matter and the gauge sector are 
effectively decoupled. 

Of course, the present study can be improved in many 
ways. Currently we are including the contributions to 
the RG flow of the four-fermion coupling arising from 
1PI diagrams with one- or two internal gluon lines [7lj . 
This will eventually allow us to get rid of the parameter 
A^ v so that we are left with a single parameter for the 
gauge and the matter sector, namely a s at a given UV 
scale. 

While a study of the QCD phase diagram in the (T, /i) 
plane is of great phenomenological importance, we think 
that a study of the phase diagram in the (T, f^) plane (or 
equivalently in the (T, | {^Tp) l 1 ^ 3 ) plane) may provide us 
with important insights concerning the interplay of chi- 
ral and confining dynamics. In particular, such a study 
of the (T, /,r) phase diagram does not suffer from prob- 
lems arising, e. g., from a complex- valued spectrum of the 
Dirac operator as it is the case at finite quark chemical 
potential \i. Therefore, an analysis of the (T, f^) phase 
diagram may be also helpful to benchmark results from 
continuum approaches against those from lattice QCD 
simulations. 
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Appendix A: Threshold functions 

In this appendix we summarize technical details con- 
cerning the derivation of the RG flow equations. 

In the computation of the RG flow equations a reg- 
ulator function needs to be specified which determines 
the regularization scheme (52j. In the present work we 
have employed an optimized spatial regulator function 
for the bosonic as well as for the fermionic degrees of 
freedom [72Tj76j . For the bosons, we choose 

Rb(p 2 )=P 2 0(k 2 -f) = A* , (Al) 

whereas we choose 

R*W=f(J^-l\ 6(k 2 -p 2 )=p ri , (J) (A2) 
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for the fermionic degrees of freedom. In the following 
we define the threshold functions relevant for the present 
work. These functions represent the 1PI diagrams con- 
tributing to the RG flow of the studied couplings. For a 
generalizations of the threshold functions to an arbitrary 
number of space-time dimensions we refer the reader to 
Ref. [Til where also the dependence of these functions on 
the anomalous dimensions is displayed. The latter is of 
no relevance for the present study. 

In order to define the threshold functions, it is con- 
venient to define dimensionless propagators for the 
bosons (B) and the fermions (tjj), respectively: 



and 



wher 



Gb(x , oj) 



G^(x ,lo) = 



1 



x + x(l + r B ) + ui 



xq + x(l + r^,) 2 + oj ' 



(A3) 



(A4) 



p 2 /k 2 



First, we define the threshold functions which appear 
in the RG flow equations for the bosonic self-interactions. 
For the purely bosonic loops, we find 



1 {t,uj) 



OO 

E 

n— — oc 



l 



dxx 2 {d t r B ) G b {u>1,uj) 



±+n B (r,w)] (A5) 



where t = T/k denotes the dimensionless temperature 
and oj n — 27rnr denotes the dimensionless bosonic Mat- 
subara frequencies. The function h B represents the Bose- 
Einstein distribution function: 



Ub{t, oj) 



1 



(A6) 



Bosonic threshold functions of order n are then derived 
from Eq. (| A5|) by taking derivatives with respect to the 
mass parameter uj: 



9 1 ( ^ 



-(n + S nt o) l n+ i(r,uj) 



(A7) 



For the purely fermionic loops contributing to the flow 
equations of the bosonic self-interactions but also to the 
RG flow of the four-fermion coupling, we find 



(F) 

l {t,u,h) 



1 1 



E 



dxx'i (d t r^)(l + ?>) x 

xG</,((j>„ + 2ttt^) 2 1 uj) 



3^1 



(l-n^,(r, i/i, w)~n^(r, —ifj,,u))) . (A8) 



Here, we have introduced the dimensionless fermionic 
Matsubara frequencies v n = (2n + 1)itt. The function 
n^j represents the Fermi-Dirac distribution function: 

1 



(A9) 



Higher-order fermionic threshold functions can again be 
found by taking derivatives with respect to the mass pa- 
rameter uj: 

^\t,u,h) = -{n + S n , )lSi(r,cJ,u). (A10) 

To prove Eq. (|20p it is convenient to perform the sum 
over Matsubara frequencies. The threshold function z£ 
can then be written as follows 

C°° 1 
4 F) (r,0, M ) = / dxxHd t r^)(l+r^\ 

Jo l4[/(x)]2 



1 



2re t [cos(27r/x)] 



e^(v^M(e 



3 I I cos(27r^) + re" 



2V/Me- 



(All) 



where f(x) = x(l+r^) 2 and 



g{x,n) = 2c 



cos(27r/i) + c 



— +1. (A12) 

The first term in the curly bracket on the right-hand side 
of Eq. (|A11|) is the zero-temperature contribution which 
is strictly positive. The second term (square bracket) 
corresponds to the finite-temperature corrections and is 
strictly negative, provided that cos(27r^) > 0. For reg- 
ulator functions with (c^r.^) > we then find that the 
finite-temperature corrections yield a negative contribu- 
tion to l{\ In the RG flow equation of the four-fermion 
coupling we sum over the eigenvalues (i = vi\4>\. Since 



E 



cos(2iwi\<f>\) > 



(A13) 



for N c = 2 and N c = 3, we have proven Eq. ([20]) . 

Finally we give the definition of the threshold function 
which appears in the RG flow equations of the Yukawa 
coupling. We have 



'i,i {t,^,h,ujb) 



E 



dx 



■*d t + 27TT/i) 2 ,W</,) x 

xG b (^,wb)}. (A14) 

To evaluate the integral over x (spatial momenta), we 
have to take derivatives with respect to the regulator 
function. For the regulator functions (|A1|) and (|A2[) these 
derivatives are given by 

l ,6{l-x)4r, (A15) 



0, 



x 1 / 2 dr^ 

-9(1- x)^-, 
x ore 



(A16) 
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where the first and the second line defines how the for- 
mal derivative d t acts on the fermion propagator and the 
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